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Abstract 

We  find  explicit  expressions  for  the  concentration  and  potential  as  functions  of  time  during  discharge  in  three  different  models  of  a  cell 
with  a  porous  metal  hydride  electrode.  The  effect  of  convection  is  studied  and  estimates  of  the  convergence  rate  to  the  steady  state  are  given. 
An  example  of  an  electric-vehicle  battery  application  is  presented. 
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1.  Introduction 

Porous  metal  hydride  electrodes  are  nowadays  widely  used 
in  electrochemical  cells.  A  battery  with  a  metal  hydride  elec¬ 
trode  has  high  energy  density  and  for  instance  in  an  electric- 
vehicle  application  one  of  the  greatest  advantages  of  metal 
hydride  batteries  is  the  low  weight  compared  to  the  traditional 
lead/acid  battery.  See  Ref.  [  1  ]  for  an  overview  of  metal 
hydride  electrodes. 

A  typical  metal  hydride  electrode  consists  of  spherical 
metal  hydride  particles  imbedded  in  a  conducting  material 
resulting  in  a  porous  structure.  When  the  electrode  is 
immersed  in  the  electrolyte  the  void  volume  is  filled.  The 
operation  of  a  metal  hydride  electrode  in  a  battery  is  based 
on  two  reactions.  The  reaction 

— Hads 

is  a  dehydridization  reaction,  where  hydrogen  is  released 
from  the  hydride  phase  forming  adsorbed  hydrogen  atoms. 
Current  is  generated  in  the  electrochemical  reaction 

xM  —  Hads+xOH~  — »xe~  +xH20-rxM 

where  the  adsorbed  hydrogen  reacts  with  hydroxide  ions  pro¬ 
ducing  electrons  and  water.  Combining  these  two  we  get  the 
total  reaction 

MH,+xOH~  ->x«-  +  xH20  +  xM 

Although  metal  hydride  alloys  and  electrodes  have  recently 
been  widely  studied  in  the  literature  there  are  only  a  few 
models  concerning  their  electrochemical  behavior.  In  Ref. 
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[2]  a  complete  cell  with  a  flat  plate  hydrogen  electrode  is 
theoretically  analyzed.  A  model  with  a  thick  metal  hydride 
electrode  is  presented  in  Refs.  [3,4]  but  this  model,  however, 
ignores  the  diffusion  and  convection  in  the  electrolyte  and 
constant  concentration  and  steady  state  throughout  the  cell  is 
assumed.  In  reality  these  assumptions  are  not  necessarily 
valid.  For  example,  the  porous  material  in  the  cell  prevents 
spontaneous  mixing  of  the  electrolyte  which  means  that  con¬ 
centration  variations  can  exist  even  in  the  steady  state.  Since 
the  concentration  is  one  of  the  crucial  factors  affecting  the 
corrosion  and  the  life  time  of  a  cell  it  is  important  to  know 
the  concentration  distribution  inside  the  electrodes.  To 
increase  the  life  time  of  electrodes  we  want  to  find  methods 
to  prevent  too  high  concentrations.  One  way  to  achieve  this 
is  to  employ  convection  through  the  cell.  Convection  can  of 
course  be  produced  externally  with  different  arrangements, 
but  a  much  more  interesting  case,  however,  is  the  natural 
convection  inside  the  cell  caused,  forinstance,  by  temperature 
gradients. 

In  this  paper  we  consider  mass  transfer  in  the  presence  of 
a  porous  metal  hydride  electrode  and  convection  parallel  to 
the  motion  of  ions  in  the  cell.  Electrode  phenomena  and 
interfacial  effects  such  as  electric  double  layers  and  surface 
overpotential  are  neglected. 

Three  different  electrochemical  systems  are  described.  In 
all  these  systems  a  porous  thick  metal  hydride  electrode  is 
placed  next  to  a  thin  counter  electrode.  The  rest  of  the  system 
is  filled  with  passive  medium  having  the  same  porosity  and 
tortuosity  as  the  metal  hydride.  The  three  cases  are;  (i)  closed 
cell  with  no  convection;  (ii)  closed  loop,  and  (iii)  an  infinite 
tube.  In  the  last  two  cases  we  include  convection. 
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In  each  case,  we  derive  explicit  expressions  for  the  con¬ 
centration  and  potential  as  functions  of  time  as  well  as  in  the 
steady  state.  The  convergence  rate  to  the  steady  state  is  esti¬ 
mated  and  the  effect  of  convection  is  investigated. 

In  a  forthcoming  article  we  include  reaction  dynamics  of 
the  metal  hydride  electrode  with  hydrogen  diffusion  in  metal 
hydride  particles  and  optimize  the  structure  of  the  porous 
electrode. 


2.  Methods  and  models 

2.1.  Preliminaries 

We  use  the  Nemst — Planck  equation  [5]  to  describe  the 
flux  density  j,  of  an  ion  species  In  the  presence  of  a  porous 
medium  with  porosity  e  the  equation  takes  the  form 

,1  5  free 

j,  =  ~e'  5D,f,cc  Vc, - — —  V<t>+  v,c,  ( 1 ) 

F 

where  c,  is  the  concentration  of  the  ion  species  i  in  the  solution 
phase,  v,  is  the  superficial  bulk  fluid  velocity  and  <t>  is  the 
potential  of  the  electric  field.  D?’QC  is  the  diffusion  coefficient 
and  <T-rco  is  the  signed  ionic  conductivity  of  the  species  i  in 
the  absence  of  the  porous  medium.  We  use  the  factor  e1  5  to 
compensate  the  effects  of  porosity  and  tortuosity  as  in  Ref. 
[6] .  The  terms  on  the  right-hand  side  of  Eq.  ( 1 )  model  the 
contributions  of  diffusion,  migration  and  convection,  respec¬ 
tively.  The  migration  term  is  often  presented  as 

^migration  =  _  ^  Ffi1  S^froc^  V(j>  (2) 

where  z,F  is  the  charge  per  mole  and  is  the  mobility  of 
the  species  i  without  the  porous  medium.  Both  forms  of  the 
migration  term  will  be  used  later. 

To  keep  notation  short  we  introduce 

D,  =  e'  5d;,cc 

or,  =  e'Vr“ 

and 

Uj  —  e1  ^uf"8 

Strictly  speaking,  the  Nernst-Planck  equations  should  be 
used  to  describe  mass  transfer  only  in  dilute  solutions  but  we 
apply  them  due  to  their  simplicity  and  the  fact  that  they  are 
widely  used  in  electrochemical  literature. 

We  consider  a  binary  electrolyte  of  two  species  of  charged 
ions  with  opposite  unit  charges  and  nonionized  solvent.  We 
make  the  customary  electroneutrality  assumption  and  denote 
the  common  concentration  by  c.  The  superscripts  ‘  +  ’  and 
‘  —  ’  of  the  diffusion  coefficients,  ionic  conductivities,  mobil¬ 
ities  and  convection  velocities  refer  to  positive  and  negative 
ions,  respectively.  From  now  on  we  use  a  convention  where 
both  ionic  conductivities  are  positive.  The  mass  balance  equa¬ 
tions  can  then  be  written  as 


ec,=  —  V •  ^  —  D  Vc  +  -^-Vd>+w  cj  +  /f  (3) 

ec,=  -  V Vc— -p-  V<£+  w+cj  +  /?+  (4) 

where  the  inhomogeneous  terms  R +  and  R~  are  production 
rates  per  unit  volume  of  positive  and  negative  ions,  respec¬ 
tively.  Generally,  the  diffusion  coefficients,  ionic  conductiv¬ 
ities,  mobilities  as  well  as  convection  velocities  depend  on 
concentration. 

Our  aim  is  to  use  analytical  methods  as  long  as  possible 
and  thus  we  simplify  Eqs.  (3)  and  (4)  considerably.  First  of 
all,  we  assume  that  the  concentration  and  electric  potential 
vary  only  in  the  axial  direction  of  the  cell.  To  keep  the  treat¬ 
ment  simple  we  use  the  velocity  of  the  solvent  as  the  convec¬ 
tion  velocity  for  all  ion  species.  This  is  a  good  approximation 
for  dilute  solutions.  Furthermore,  we  ignore  the  viscous 
forces  and  assume  that  the  solvent  velocity  is  axial  and  con¬ 
stant  throughout  the  cell.  The  (scalar)  solvent  velocity  is 
denoted  by  v.  Hence  Eqs.  ( l)-(4)  reduce  to  a  one-dimen¬ 
sional  form.  By  assuming  that  the  porosity  is  constant 
throughout  the  cell  we  avoid  discontinuities  in  the  coefficients 
of  our  mass  balance  equations.  The  effect  of  axial  dispersion 
on  diffusion  coefficients  has  been  neglected.  Finally,  we 
assume  that  the  cell  is  discharged  with  constant  current  and 
that  all  the  coefficients  and  the  initial  state  are  constant. 

We  apply  two  different  models  for  the  migration  term. 
When  the  concentration  variations  in  the  cell  are  small  or  a 
simple  model  is  desired  we  use 

j- migration  _  (5) 

where  the  ionic  conductivities  tr  *  are  constants.  If  the  con¬ 
centration  is  expected  to  vary  significantly  it  is  better  to  define 

^mtgrauon  _  +FU  ±C<}>X  (6) 

where  the  mobilities  are  constant  but  there  is  a  linear 
dependence  on  concentration. 

When  modelling  electrochemical  cells  it  is  customary  to 
assume  that  there  are  concentration  variations  only  in  thin 
boundary  layers  next  to  the  electrodes  and  that  the  bulk  of 
the  electrolyte  remains  in  the  constant  initial  state.  In  our 
application,  however,  porous  medium  fills  the  whole  cell  and 
thus  we  assume  that  there  are  neither  spontaneous  mixing  nor 
other  processes  that  would  even  out  concentration  differences 
in  the  electrolyte.  Hence  we  apply  Eqs.  (3)  and  (4)  through¬ 
out  the  cell,  not  only  near  the  electrodes. 

Usually  the  inhomogeneous  terms  R*  and  R~  arise  from 
chemical  reactions  taking  place  in  the  electrolyte.  Now  there 
are  no  chemical  reactions  in  the  electrolyte  itself  but  OH~ 
ions  are  produced  at  the  counter  electrode  and  destroyed  on 
the  surface  of  the  porous  metal  hydride  electrode.  Thus  we 
can  use  R  —  /?sourcc  f^smk  2nd  R  ^source  f^sink  1®  model 
the  electrodes  as  in  Ref.  [6] . 

Accordingly  the  porous  metal  hydride  electrode  appears 
as  a  sink  term  for  OH  ~  ions  in  Eq.  ( 3 )  in  all  our  three  cases. 
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On  the  other  hand,  it  is  less  conventional  to  use  the  same 
approach  in  describing  the  thin  counter  electrode  with  a 
source  term  instead  of  a  boundary  condition.  This  approach 
is  very  attractive  due  to  its  simplicity.  We  do  not  have  to  pay 
attention  to  directing  the  ions  forming  at  the  counter  electrode 
to  up  or  down  stream.  Instead,  the  ion  flux  from  the  counter 
electrode  is  controlled  by  diffusion,  migration  and  convection 
in  the  same  way  as  everywhere  else  in  the  cell. 

Since  the  positive  ions  are  passive  and  do  not  take  part  in 
reactions  we  have  R+  =0  everywhere. 

During  discharge  hydrogen  is  consumed  according  to  the 
local  reaction  rate  which  varies  throughout  the  electrode.  The 
reaction  rate  depends,  among  other  things,  on  the  concentra¬ 
tion  and  potential  of  the  electrolyte,  the  potential  of  the  metal 
hydride  and  the  concentration  of  the  hydrogen  atoms  on  the 
surface  of  the  metal  hydride  particles  [4].  In  optimizing  the 
performance  of  the  electrode  one  goal  is  to  ensure  that  hydro¬ 
gen  is  consumed  uniformly.  Otherwise  some  parts  of  the 
electrode  would  become  passive  before  the  total  hydrogen 
storage  is  exhausted  and  thus  the  efficiency  of  the  battery 
would  decrease  towards  the  end  of  the  discharge. 

Thus  we  assume  that  all  the  OH  ~  ions  that  are  produced 
on  the  counter  electrode  react  on  the  metal  hydride  electrode 
and  that  this  reaction  rate  is  uniform  throughout  the  metal 
hydride,  corresponding  to  the  optimal  operation  of  the 
electrode. 

Hence,  if  the  counter  electrode  is  located  at  x= jc,  and  the 
metal  hydride  electrode  occupies  the  region  between  the 
points  x = X|  and  jc  =  x, ,  we  set 


if  jc^O 
if  x>0 


and  abbreviate /»,,(*)  =h(x—y). 

Note  that  according  to  our  assumptions  the  ratio  of  the 
ionic  conductivities  of  positive  and  negative  ions  is  constant. 
Then  the  potential  gradient  can  be  eliminated  from  Eqs.  (3) 
and  (4)  and  the  concentration  depends  only  on  the  ratio  of 
the  migration  terms  but  not  on  their  explicit  form. 

In  all  models  presented  below  the  negative  ions  move  from 
left  to  right  in  the  metal  hydride  electrode  so  >'\at  the  current 
is  in  the  opposite  direction. 


2.2.  Closed  cell  without  convection 


We  consider  a  cell  with  a  porous  metal  hydride  electrode 
filling  the  region  (0,  /]  and  a  thin  counter  electrode  placed  at 
jc = 0,  see  Fig.  1 .  We  assume  that  there  is  no  flow  through  the 
boundaries  at  jc= 0  and  x— l  and  that  there  is  no  convection 
in  the  cell  so  that  only  diffusion  and  migration  contribute  to 
the  mass  transport. 

Then  the  concentration  and  potential  satisfy  the  equations 


r_ 

/ 


(7) 

(8) 


for  0  <  jc  <  1,  t  >  0  with  the  boundary  conditions 


/(sink  "  J 


—  JTcl 


where 


Xd(x)  = 


if  x&D 
otherwise 


is  the  characteristic  function  of  the  interval  D  and  J~  is  the 
current  density  at  the  counter  electrode  expressed  in  flux 
density  units  mol/m2  s.  Then 


J  *ik(0 


(-d  '.+v*)L„--r 

(-dv--T*)L-o 

(-°‘c'+T*)L,-o 

that  hold  for  all  t>  0.  The  initial  condition  is 
c(  0,  x)=c0 


(9) 

(10) 

(11) 

(12) 


(13) 


as  desired.  In  the  closed  cell  case  the  counter  electrode  is 
modelled  as  a  boundary  condition 

so  that  /(SoUrcc  =  0  but  in  the  other  two  cases  we  set 

^source-  7  ^x\ 

where  6fl(x)  =  8(x  —  jc,)  is  a  Dirac  delta  functional  centered 
at  x = jq.  That  is,  the  counter  electrode  is  described  as  a  point¬ 
like  source  of  OH-  ions. 

In  the  following  sections  we  frequently  use  the  Heaviside 
step  function  h  defined  by 


forO  <,x^l. 


Fig.  1 .  Closed  cell. 
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Note  that  according  to  the  boundary  conditions  (9)  and 
( i  i )  and  the  inhomogeneous  sink  term  in  Eq.  (7)  the  current 
fed  into  the  electrolyte  at  x  =  0  is  uniformly  gathered  by  the 
porous  metal  hydride  electrode  on  the  interval  (0,  /)  so  that 
the  flux  of  the  negative  ions  vanishes  at  x=/.  On  the  other 
hand,  the  boundary  conditions  (10)  and  (12)  and  Eq.  (8) 
are  formulated  so  that  the  positive  ions  do  not  react  on  the 
electrodes. 

For  later  use,  we  derive  an  equivalent  set  of  boundary 
conditions  by  solving  the  Eqs.  (9)-(  12)  for  the  concentra¬ 
tion  and  potential  gradients  at  the  boundaries: 


cr+J- 


D'cr~  +  D~  cr^ 


cx\x~,  =  0 


FD+J~ 


Mx-a  D  +  gr-+D-(T  + 

<kL-;  =  0 


(14) 

(15) 

(16) 
(17) 


From  Eqs. (7)  and  (8)  and  the  boundary  conditions  (14) 
and  (15)  we  obtain  after  some  manipulations 


ec,=DCxx— 


t+r 


with 


Cjc'j[  =  0~ 


t+J~ 


o 

where 

„  D~cr*  +D  +  cr~ 
cr+  +a~ 

is  the  integral  diffusion  coefficient  and 


(18) 


(19) 


(20) 


t  = 


tr+  +<r 


is  the  transference  number  of  positive  ions. 

In  order  to  use  the  separation  of  variables  method  [7]  for 
solving  the  system  Eqs.  ( 1 8)— (20)  with  the  initial  condition 
(13)  we  make  a  change  of  the  unknown.  We  define 

t*J-  ,  t+J- 
u-c-  — —  jtH — —x 
2  ID  D 

Then  the  original  problem  transforms  to  the  homogeneous 
diffusion  equation 

eu,=Ou„  (21) 

with  the  homogeneous  bouno.ry  conditions 

«,L-o=«xl*-i=0  (22) 

and  the  initial  condition 


t*j-  ^  t+J~ 

w(0,  x)  =  c0  - _ ■  xr  3 — —x 


2  ID 


D 


(23) 


The  transformed  problem  (21 )— (23)  is  easily  solved  and  the 
solution  of  the  original  problem  is  then  given  by 

f+/-/ri  x  lx2 
c«.*)-c„+—  [5-7+5?- 

“1  /  n2ir2D\  (mt \~l 

2L^exp(— f?t'MtxJJ  (24) 

Here,  as  well  as  in  the  two  other  cases,  the  partial  differential 
equation  can  alternatively  be  solved  with  Laplace  transforms 
but  the  solution  is  again  expressed  in  terms  of  a  series. 

With  the  constant  ionic  conductivities  (5),  by  subtracting 
Eq.  (8)  from  Eq.  (7)  and  rearranging  the  terms,  we  obtain 


_a 

3x! 


(D--D+)Cxx—j- 


(25) 


Then,  integrating  Eq.  (25)  twice  with  respect  to  x  and  using 
the  boundary  conditions  ( 14)— ( 17)  yields 


<#>(',  x)=— -7— r[ 

cr  +<r+  L 


(D~  —  D+ )c(t,  x) +J 


H] 


(26) 

Similarly,  using  the  linear  concentration  dependence  (6) 
for  the  ionic  conductivities  leads  to 


c(t,x)_ 
0) 


(27) 


From  Eq.  (24)  it  is  easy  to  see  that  the  steady  state  solution 
defined  by  C(x)  =lim,_„c(t,x)  is 


s  t+J~l(l  x  lx*) 

CM-c°+—(n*2P) 

and  that  for  all  xe  [0,/] 
lc(,.<)  -C(X)| 


(28) 


4  t+J-l  ( 
=~^exp {' 


ipD  \ 

7? 


That  is,  the  concentration  converges  uniformly  to  the  steady 
state  with  exponential  rate.  Note  that  shortening  the  electrode, 
i.e.  making  l  smaller,  makes  the  convergence  faster. 

Expressions  for  the  steady  state  potentials  can  be  derived 
by  substituting  the  steady  state  concentration  given  by  Eq. 
(28)  into  Eqs.  (26)  and  (27). 

2.3.  Circular  cell  with  convection 

The  reason  for  studying  a  circular  cell  is  that  we  wanted  to 
find  a  simple  model  compatible  with  forced  or  natural  con¬ 
vection  through  the  electrodes.  The  opposite  ends  of  a  closed 
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Porous  metal  hydi&e  electrode 


Fig.  2.  Circular  cell. 


cell  like  the  one  described  previously  are  connected  with  a 
tube  making  electrolyte  circulation  through  the  electrodes 
possible.  Below  we  present  an  idealization  of  this 
arrangement. 

We  consider  an  imaginary  electrochemical  system  con¬ 
sisting  of  a  thin  tube  in  which  there  are  two  electrodes 
(Fig.  2).  The  metal  hydride  electrode  occupies  the  interval 
(jt„  jc2]  and  the  infinitely  thin  counter  electrode  is  placed  at 
x=x,.  The  ends  of  the  tube  are  connected  and  it  forms  a 
closed  loop.  Thus  there  are  two  routes  for  the  current  to  pass 
from  one  electrode  to  the  other. 

The  mathematical  model  is  as  follows.  The  concentration 
c  and  potential  <t>  satisfy  the  equations 


Finally  we  set  an  initial  condition 

c(0,r)=c0,  0  £x<L  (35) 

Eliminating  the  migration  terms  from  Eqs.  (29)  and  (30) 
leads  to  an  equation  for  the  concentration 

(36) 

\  X2-X,J 

The  periodic  boundary  conditions  (31)  and  (34)  encour¬ 
age  us  to  use  a  Fourier  series  method  (see  Ref.  [7] )  to  solve 
the  problem  involving  Eqs.(36),  (31),  (34)  and  (35).  We 
first  solve  the  problem  with  zero  initial  condition  and  then 
add  c0  to  the  solution  to  find  the  solution  to  the  original 
problem. 

We  define 

C„(0=  £  J  c(  r,  £)  exp^  -  i^pfj  df 
0 

where  n  is  an  integer.  Then  Eq.  (36)  transforms  to 

_ ,  l  An2 -n2  .2n7r\ 

eC'Jt )  =  -(d— ^+uri— jCn(r) 

+t+r{r;exp(-i?fJCi) 

“2^-*,) 


ec'=-s("D+c'~V<k+cw) 


(29) 

(30) 


forOCjc^Land  f>0. 

Due  to  the  circular  nature  of  the  cell  the  boundary  condi¬ 
tions  are  periodic: 

c(r,0)=c(f,L)  (31) 

<fe(/,0)  =  <k(f,  L)  (32) 

and 

j-(t,0)=j-(.t,L),  j+{t,Q)=j+(t,L)  (33) 

where 


j~  —  —D  cx+—<f>x+cv 


whose  solution  with  the  transformed  initial  condition 
C„(0)=0is 

„  +  r An2™1  ,2nii\  '  f  1  J  .2 nv  \ 

c„(o=t  j  J  \zexpV-,~ c,j 

i  r  (  ,2mr  \  (  .2nir  'll! 

-2***-*.)  H_,Txr 

{'-eXP[-(^ 


“^77^  ,2vMt\  t 

-+i — : 


for  n  0. 

We  r.ow  determine  C0.  It  follows  from  Eq.  (30)  and  the 
boundary  condition  (33)  that  the  mass  is  conserved 

L 


>■ 


0df=O 


With  the  temporary  zero  initial  condition  this  means  that 


and 

j+  =  ~D+cx-~-<l>x  +  cv 

In  other  words,  we  require  that  the  concentration,  electric 
field  and  flux  density  are  periodic.  It  follows  from  Eqs.(31  )— 
(33)  that  the  concentration  derivative  is  periodic  as  well 

(34) 


L 

Jc(r,£)d£=0 

0 

for  t  S  0.  Thus  C0(  t )  =  0. 

Then  the  solution  to  Eq.  (36)  together  with  the  boundary 
and  initial  conditions  is  given  by  a  doubly  infinite  series 

ao 

c(r,x)=c0+  £  C„(t)  ex 

n“  —  00 


cx(t,  0)  =  c,(f,  L) 
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Since  the  concentration  is  a  real  function  we  can  also  express 
it  as: 


c(t,x) 


=  c„+  £  an 

71  =  |L 
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(2mr  \ 

IT'J 


+  b„(t)  sin 
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(39) 

To  find  an  expression  for  the  potential  when  the  ionic 
conductivities  are  constant  (model  (5) )  we  subtract  Eq.  (29) 
from  Eq.  (30)  to  eliminate  the  time  derivatives  and  convec¬ 
tion  terms  and  solve  for  <£„.  This  yields 

7™r  [  (O' 

<T  +<T  L  V  X2-X,/J 

Integrating  the  previous  equation  with  respect  to  x  leads  to 

4>At, X)  — :  T (D-  — D+)cx(t . X) 
cr  +cr  L 

X2  X,  J 

where  A  is  an  arbitrary  function  of  t.  Note  that  this  expression 
is  periodic  for  every  choice  of  the  integration  constants.  To 
determine  the  potential  we  must  thus  find  an  extra  condition 
for  <f>.  It  is  very  natural  to  require  that  that  the  potential  itself 
is  periodic  too 
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for  all  f>0.  With  this  choice  we  get,  after  an  additional 
integration. 
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(40) 


When  the  ionic  conductivities  depend  linearly  on  concen¬ 
tration  (model  (6))  the  potential  can  be  expressed  as 
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4>(i,  x) 


F(u+  +u~) 


j(Z>~ -D+)\nc(t,x) 
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f _ hzi _ 

Jc(t,0(x2-X,) 

0 


where  If  is  a  constant.  That  is,  the  solution  converges  uni¬ 
formly  to  the  steady  state  with  exponential  rate.  As  before, 
decreasing  the  elect- :>lyte  volume,  i.e.  making  L  smaller, 
increases  the  convergence  rate. 

Expressions  for  the  steady  state  potentials  can  be  derived 
by  substituting  the  steady  state  concentration  given  by  Eq. 
(42)  intoEqs.  (40)  and  (41). 
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From  the  series  representation  (37)-(39)  it  is  easy  to  see 
that  there  is  a  steady  state  solution  C(x)  =  lim, -,»c(r,  x) 
given  by 
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2.4.  Infinite  tube  with  convection 

This  is  a  limiting  case  of  the  circular  cell.  If  we  let  the 
length  of  the  loop  tend  to  infinity  we  obtain  a  straight  infinite 
tube  as  the  limit.  This  case  is  of  interest  because  it  can  be 
thought  to  model  a  cell  connected  to  two  big  reservoirs  of 
electrolyte.  Contrary  to  the  circular  cell,  there  is  no  secondary 
route  for  the  current.  That  is,  all  the  current  passes  through 
the  interval  between  the  electrodes. 

An  infinite  tube  is  filled  with  electrolyte,  a  porous  metal 
hydride  electrode  occupies  the  region  ( 0,  /]  and  a  thin  counter 
electrode  is  placed  at  x=0  (Fig.  3).  The  concentration  and 
potential  satisfy  the  equations 

ec,=  ^-D'c,+— ^  +  ucj+y^fiD-^yilj  (43) 

ec,=  -^(-Z>+c*-^r4>,  +  rc)  (44) 

for  -  »  <  x  <  <»,  t  >  0  with  the  initial  condition: 

c(0,  x)  -c0  (45) 

for  —  co<x<°°. 

We  use  Fourier  transformations  (see  Ref.  [7])  to  solve 
the  problem  involving  Eqs.  (43)--(45).  Eliminating  the 
migration  terms  from  Eqs.  (43)  and  (44)  results  in 
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and 

b„(°°)  =  lim  b„(t) 


f  1  2nir  1 

+  vL  -  cos-— x,  -  - — - - - 

Li  i  2nir(x2-X|) 

(  .  2nir  .  2 ntr  YT\ 
x(sm— r2-sm—x,JJ| 

[0.1]: 

(  4 Dtr2 


Furthermore,  for  all  re  [0,  L] : 
\c(t,x)-C(x)  \  <Kex\ 


where  D  and  t+  are  lefined  as  before.  The  solution  of  Eq. 
(46)  with  Eq.  (45)  is  given  by 


t+J~L2  j 

r2Dn\ 

1  .  2n7r 

c(r,  x)  —  c0  = 

4D2n3  7T3  + 1 rL2n  -tr  1 

2Dn  ,t| 

jsm~rx> 

1  l 

2nir 

2mr  \~l 

J  [x 

expS - 

+  1  CO! 

2nir(x2-xl)  \ 

nr* 

-cos— X.JJ 

r  1 

t/e  w 

-c0=  I"  f  rr 

J.  J  j4irD(t/e—  r) 


4D(t/e—  r) 


-}[ 


«b(f) 


t/f 

-jwejdf*. 

L  4 Dt  J  2/|_  V  y4Dr  / 

Ms)]}- 

where  erf(  • )  is  the  error  function. 


(47) 


34 


J.  Heikonen  et  al.  /Journal  of  Power  Sources  62  (1996)  27-39 


Paash/e  peraua  material 


Fig.  3.  Infinite  tube. 

Let  us  now  take  a  look  at  the  potential  when  the  ionic 
conductivities  are  constant  (model  (5)).  We  subtract  Eq. 
(43)  from  Eq.  (44)  and  solve  for  This  yields 

We  integrate  twice  with  respect  to  x  and  require  that  the 
electric  field  vanishes  at  inifinity.  As  the  result  we  get 

<Ht ,  x)  =—  £-  —  { (ZT  - D+  )c(t,  x) 
a  +cr  (. 

+ ^  -  n  (*> + ih,M  ]}  (48) 

Suppose  now  that  the  ionic  conductivity  depends  linearly 
on  concentration  (model  (6)).  Then  the  potential  is  given 
by 

0  0 

(49) 

Next  we  investigate  the  steady  state.  We  let  C=lim,_„c 
and  define 

f~  Sq— yAfto. /i 


so  that 


ey  i-  f .,_/(l-Jt)//  if  .re  (0,  /] 
\  0  if  jcS  [0,  /] 

Then,  using  partial  integration,  we  find  that 
C(x)  -c(l,x) 

f  1  [  (x-£-vr)2~ 

l  J^expL 

00 


4Dr 

(x-g-vr)2 

4Dt 


/(g)  dgdr 


]f(g)i|:". 


f2((J c-S-w)  r  (*-g-vr)2l_  ,,1  . 

-  J - adt - expL - J F(f)  df; dT 

0 

_ , +r  f  1  f2((x-g)-vr) 

■j  tJAttDt)  4Z>t 

z/e  0 

* exp[ ~ (* ] F( A d£ ' dT 


x{Sign(rt-™r[(J,“^"*]}d{ 


From  this  it  is  easy  to  see  that  the  solution  converges  point- 
wise  to  the  steady  state.  Now  assume  that  x  e  [0,  1] .  If  v>  0, 
then 

ICM 

For  t/<  Owe  have 

Finally,  if  v=  0  ii  follows  that 

|C(x)-c(r,x)l5^-/erf(^=) 

Thus,  the  larger  the  value  of  I  v|  the  faster  the  rate  of 
convergence. 

It  is  clear  from  Eq.  (47)  that  the  steady  state  solution  is 
given  by 

M^sr)  ■  "(is;)]}  dr  <501 

Expressions  for  the  steady  state  potentials  can  be  derived 
by  substituting  the  steady  state  concentration  given  by  Eq. 
(50)  intoEqs.  (48)  and  (49). 


3.  An  application 

3.1.  Parameters 

In  all  our  examples  KOH  is  used  as  the  electrolyte  solution. 
The  initial  concentration  is  chosen  to  be  5  M.  The  ionic 
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Table  1 

Electrochemical  data 


Parameter 

Value 

Refs. 

Dfree  (m2/s) 

3.23X10-’ 

[8] 

4"“(m2/flequiv) 

1.17X10-2 

[8] 

t  + 

0.22 

[9] 

f/r  (m2/s) 

2.07X10-’ 

If™  (m2/s) 

7.34X10-’ 

e'rd/nm) 

12.87 

^(l/flm) 

45.63 

(m2mol/J  s) 

2.76X10-’3 

u?*  (m2mol/J  s) 

9.80X10'13 

F  (C/equiv) 

96500 

diffusion  coefficients,  conductivities  and  mobilities  for  5  M 
KOH  electrolyte  are  derived  from  the  integral  diffusion  coef¬ 
ficient  and  the  equivalent  conductance  using  the  transference 
numbers.  The  values  for  these  parameters  are  given  in 
Table  1. 

The  absolute  value  of  the  convection  velocity  |  ul  varies 
between  0  and  10“  5  m/s.  The  velocities  used  here  can  easily 
be  shown  to  be  reasonable  and  are  in  good  consistency  with 
even  quite  small  temperature  differences,  10  °C  for  instance, 
frequently  occurring  in  electrochemical  systems. 

Next  we  discuss  the  parameters  describing  the  cell,  the 
porous  electrode  and  operating  conditions  that  are  given  in 
Table  2. 

The  cell  parameters  are  chosen  for  electric-vehicle  appli¬ 
cations,  for  example.  Then  it  is  convenient  to  assume  a  dis¬ 
charge  time  of  about  3  h  and  a  discharge  current  density  of 
100  mA/cm2  corresponding  to  a  flux  density  of  0.01  mol/ 
m2  s  of  OH"  ions  and  a  capacity  of  300  mAh/cm2  for  the 
electrode.  For  instance,  with  a  much-studied  metal  hydride 
alloy  MmNi3.sCoo.8Alo.4Mn0.3  [10,11]  having  the  capacity 
of  about  250-270  mAh/g  this  means  that  the  surface  density 
of  the  metal  hydride  alloy  is  1.1-1 .2  g/cm2. 

The  thickness  of  the  electrode  and  its  porosity  are  related. 
It  seems  that  the  optimal  value  for  the  porosity  is  about  0.5. 
Then,  if  the  electrode  has  3  wt.%  polytetrafluoroethylene 
(PTFE)  and  1  wt.%  carbon  as  the  electrically  conducting 
binder  material  it  follows  that  the  electrode  should  be  3  mm 
thick.  The  thickness  of  the  electrode  is  thus  relatively  high 
but  it  has  been  chosen  keeping  in  mind  sintering  fabrication 
methods  with  which  high  discharge  rates  even  for  thick  elec¬ 
trodes  can  be  achieved. 

Table  2 

Cell  parameters 


Parameter  Value 


J  (mol/m2s) 

0.01 

1  (m) 

0.003 

Mm) 

0.1 

€ 

0.5 

JC,  (m) 

0.001 

*!  (m) 

0.004 

c0  (mol/1) 

5 

With  the  thickness  of  3  mm  the  capacity  of  the  electrode 
is  1  Ah/cm3.  To  balance  this  we  can  use  an  NiOOH  electrode 
[12]  or  an  oxygen  electrode  as  the  counter  electrode,  for 
example. 

In  the  case  of  the  circular  cell  the  total  length  of  the  tube, 
L,  is  chosen  to  be  large  compared  with  the  length  of  the  metal 
hydride  electrode,  /= x2  -  xt.  In  this  way  the  secondary  cur¬ 
rent  outside  the  metal  hydride  electrode  is  negligible. 


4.  Results 

Applying  the  previously  introduced  mathematical  models 
with  the  parameter  values  given  above  we  computed  concen¬ 
tration  and  potential  profiles  in  various  situations.  We  empha¬ 
size  that  by  potential  we  mean  the  potential  in  the  electrolyte 
solution.  Also,  when  using  the  phrases  constant  and  linear 
ionic  conductivities  we  refer  to  models  (5)  and  (6),  respec¬ 
tively.  In  all  cases  below,  the  concentration  and  potential 
profiles  only  within  the  porous  electrode  are  shown,  that  is, 
we  are  not  concerned  with  the  concentration  and  potential 
outside  the  metal  hydride  electrode. 

In  Fig.  4  the  time  dependence  of  concentration  with  the 
closed  cell  model  is  presented.  After  the  discharge  time  of 
3000  s  no  change  in  the  concentration  profile  can  be  observed 
so  that  an  approximate  steady  state  has  been  reached.  In  this 
section  we  use  the  phrase  steady  state  in  this  sense.  The 
difference  between  the  maximum  and  minimum  concentra¬ 
tions  inside  the  metal  hydride  electrode  increases  with  dis¬ 
charge  time.  At  the  steady  state  a  difference  of  2.8  M  has 
developed.  The  corresponding  potential  distributions  with 
constant  and  linear  ionic  conductivities  are  shown  in  Figs.  5 
and  6,  respectively.  The  potential  difference  between  the  ends 
of  the  metal  hydride  electrode  decreases  with  discharge  time. 
At  the  steady  state  the  difference  is  45  mV  for  the  constant 
ionic  conductivity  model  and  42  mV  for  the  linear  ionic 
conductivity  model. 

To  examine  the  possibility  to  change  the  concentration 
distribution  and  especially  concentration  differences  inside 


Fig.  4.  Concentration  in  the  closed  cell  as  a  function  of  a  time. 
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Fig.  5.  Potential  in  the  cell  as  a  function  of  time  with  constant  ionic 
conductivities. 


Fig.  6.  Potential  in  the  closed  cell  as  a  function  of  time  with  linear  ionic 
conductivities. 


Fig.  7.  Concentration  in  the  circular  cell  as  a  function  of  time  with  no 
convection. 

the  electrode  by  using  convection  the  two  other  cell  models, 
i.e.  circular  cell  and  infinite  tube,  were  studied.  The  physical 
validity  and  correspondence  of  these  cases  compared  to  the 


closed  cell  were  first  investigated  by  calculating  the  concen¬ 
tration  and  potential  distributions  without  convection.  The 
concentration  distributions  for  the  circular  cell  and  the  infinite 
tube  are  presented  in  Figs.  7  and  8,  respectively.  Comparing 
these  to  Fig.  4  we  see  that  the  shapes  of  the  concentration 
distributions  are  practically  equal  in  each  case  and  that  the 
difference  between  the  maximum  and  minimum  concentra¬ 
tions  has  remained.  On  the  other  hand,  the  maxima  and  min¬ 
ima  of  the  concentrations  have  become  lower.  This  is  due  to 
the  effect  of  the  electrolyte  reservoirs  in  the  models.  The  main 
observation  is  however  that  the  time  needed  to  achieve  a 
steady  state  is  in  both  the  circular  cell  and  the  infinite  tube 
case  beyond  the  assumed  discharge  time  of  10  000  s  and  the 
times  are  now  of  the  order  of  5  X  105  and  5  X  106  s,  respec¬ 
tively.  This  is  an  obvious  consequence  of  the  larger  electro¬ 
lyte  volumes.  The  potential  distributions  in  the  circular  ceil 
and  the  infinite  tube  cases  are  not  included  because  they  are 
naturally  almost  equal  to  those  of  the  closed  cell  case. 

The  effect  of  convection  on  the  concentration  distributions 
at  the  steady  state  is  shown  in  Figs.  9  and  10  with  different 
convection  velocities.  With  low  convection  velocities  of 
±  10“7  m/s  the  differences  between  the  maximum  and  min¬ 
imum  concentrations  arc  practically  the  same  as  in  the  cases 
without  convection,  but  the  maxima  and  minima  of  the  con¬ 
centrations  have  become  either  lower  or  higher  than  without 
convection,  depending  on  the  direction  of  convection.  With 
convection  velocities  of  ±  10“ s  m/s  the  shapes  of  the  con¬ 
centration  profiles  and  differences  between  maximum  and 
minimum  concentrations  are  already  markedly  changed. 
When  the  convection  is  against  the  flow  of  negative  ions  the 
shape  of  concentration  in  the  electrode  is  still  continuously 
decreasing,  but  with  the  convection  in  the  other  direction  the 
shape  has  been  changed  much  more  remarkably,  that  is,  the 
minimum  of  the  concentration  has  moved  towards  the  counter 
electrode.  In  either  case  the  difference  between  the  maximum 
and  minimum  concentrations  has  decreased  to  1  M  from  2.8 
M  for  the  case  without  convection.  With  quite  high  convec¬ 
tion  velocities  of  ±  10“  5  m/s  the  concentration  profiles  are 


Fig.  8.  Concentration  in  the  infinite  tube  as  a  function  o»  time  with  no 
convection. 
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Fig.  9.  Steady  state  concentration  in  the  circular  cell  for  different  convection 
velocities. 


Fig.  1 1.  Steady  state  potential  in  the  infinite  tube  for  different  convection 
velocities  with  constant  ionic  conductivities. 


velocities. 

further  flattened,  so  that  the  difference  between  the  maxima 
and  minima  is  only  0.2  M.  Also  the  changes  in  the  shape  of 
the  concentration  profiles  observed  with  the  velocities  of 
±  10~6  m/s  are  further  strengthened. 

Note  that  in  Fig.  10  the  concentration  at  the  down  stream 
end  of  the  metal  hydride  electrode  in  the  infinite  tube  model 
is  equal  to  S  M,  the  initial  concentration,  when  there  is  con¬ 
vection.  This  is  due  to  the  requirement  of  boundedness  of 
solutions  of  the  mass  balance  equations.  In  the  circular  cell 
case  (Fig.  9)  the  down  stream  end  concentration  is  pretty 
close  to  S  M. 

The  potential  distributions  at  the  steady  state  with  different 
convection  velocities  in  the  cases  of  constant  and  linear  ionic 
conductivities  are  presented  in  Figs.  1 1  and  1 2  for  the  infinite 
tube  model.  The  results  for  the  circular  cell  are  not  shown 
because  they  are  practically  equal  to  those  of  the  infinite  tube. 

In  Fig.  1 1,  where  the  constant  ionic  conductivity  model 
was  used,  the  potential  distributions  corresponding  to  low 
convection  velocities  of  ±  10-7  m/s  are  very  close  to  the 
potential  distribution  without  convection.  On  the  other  hand. 


with  the  linear  model  for  ionic  conductivities,  shown  in  Fig. 
12,  these  velocities  cause  noticeable  differences.  The  poten¬ 
tial  difference  between  the  electrode  ends  changes  from  46 
mV  (no  convection)  to  37  mV  (convection  in  the  direction 
of  negative  ions)  and  to  71  mV  (convection  against  the 
negative  ions).  This  difference  between  Figs.  11  and  12  is 
easily  explained  by  the  concentration  changes  for  velocities 
0  and  ±  10~7  m/s  shown  in  Fig.  10  together  with  the  linear 
ionic  conductivity  model. 

From  Fig.  1 1  it  can  be  seen  that  the  potential  difference 
between  the  ends  of  the  metal  hydride  electrode  increases 
from  45  mV  (no  convection)  to  58  or  to  65  mV  with  the 
convection  velocity  of  10-6  m/s  along  with  or  against  the 
negative  ions,  respectively,  when  constant  ionic  conductivi¬ 
ties  are  assumed.  With  the  linear  ionic  conductivities  in  Fig. 
12  the  potential  difference  between  the  ends  of  the  metal 
hydride  electrode  has  changed  from  46  to  50  mV  or  79  mV 
when  the  convection  is  10~6  m/s  along  with  or  against  the 
negative  ions,  respectively.  Again  the  direction  of  the  con¬ 
vection  affects  the  potential  more  when  the  linear  ionic  con¬ 
ductivity  model  is  used. 


Fig.  12.  Steady  state  potential  in  the  infinite  tube  for  different  convection 
velocities  with  linear  ionic  conductivities. 
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With  velocities  of  ±  10-5  m/s  the  effect  of  the  direction 
of  the  convection  is  very  much  decreased  especially  in  the 
case  of  constant  ionic  conductivities  as  can  be  seen  from  Fig. 
1 1 .  The  potential  difference  between  the  electrode  ends  has 
increased  from  45  to  about  69  mV  regardless  of  the  direction 
of  convection.  With  the  linear  ionic  conductivity  model  in 
Fig.  1 2  the  potential  difference  has  increased  from  46  to  about 
70  mV. 


5.  Conclusions 

From  the  above  results  it  is  obvious  that  the  concentration 
distributions  cannot  be  ignored  in  thick  porous  electrodes. 
Our  examples  show  that  quite  remarkable  concentration  gra¬ 
dients  are  possible  in  such  electrodes.  It  was  also  observed 
that  the  time  needed  to  reach  the  steady  state  is  much  longer 
than  what  was  expected.  When  electrolyte  reservoirs  were 
present,  as  in  the  models  of  circular  cell  and  infinite  tube  in 
this  presentation,  this  time  was  beyond  the  supposed  dis¬ 
charge  time  of  3  h. 

The  concentration  gradients  can  be  affected  by  convection, 
if  a  suitable  velocity  is  chosen.  With  high  enough  velocities 
(of  the  order  of  10~5  m/s  in  our  examples)  the  concentration 
profiles  can  be  flattened  almost  totally,  while  with  low  veloc¬ 
ities  (of  the  order  of  10~7  m/s  in  our  examples)  no  effect 
on  the  difference  between  the  maximum  and  minimum  of 
concentrations  is  seen. 

Even  though  flat  concentration  profiles  can  be  thought  to 
be  advantageous  because  of  lower  concentration  maxima  and 
therefore  reduced  corrosion,  flat  profiles  at  the  same  time 
generally  correspond  to  higher  potential  gradients  inside  the 
electrode.  Since  the  electrolyte  potential  has  a  clear,  and  in 
some  circumstances  even  very  crucial,  effect  on  corrosion 
rate,  the  steep  potential  profiles  are  also  undesirable.  Thus  it 
seems  that  in  trying  to  extend  the  life  time  of  a  cell  optimizing 
the  concentration  and  potential  profiles  is  of  great  importance. 

However,  more  attention  should  be  paid  to  modelling  the 
concentration  dependence  of  the  ionic  conductivities.  In  this 
presentation  only  constant  and  linear  models  for  the  concen¬ 
tration  dependence  of  ionic  conductivities  were  used,  because 
of  the  requirement  of  simplicity.  Using  these  models,  quite 
different  results  were  obtained  with  lower  velocities.  With 
higher  velocities  the  models  gave  almost  identical  results 
because  the  concentration  differences  in  the  electrode  were 
very  small. 


6.  List  of  symbols 

c  concentration  of  negative  and  positive  ions 

C  steady  state  concentration  of  negative  and 

positive  ions 

c0  initial  concentration  of  negative  and  positive  ions 

D  effective  integral  diffusion  coefficient 


DUcc  integral  diffusion  coefficient 

D+  effective  diffusion  coefficient  of  positive  ions 
D  ~  effective  diffusion  coefficient  of  negative  ions 
D *5“  diffusion  coefficient  of  positive  ions 

£frec  diffusion  coefficient  of  negative  ions 

erf  error  function 

F  Faraday’s  constant 

h(  Heaviside  step  function  with  the  step  at  x  =  f 

j+  scalar  flux  density  of  positive  ions 

j ~  scalar  flux  density  of  negative  ions 

J~  discharge  current  in  flux  density  units 
l  length  of  the  porous  electrode 

L  length  of  the  cell 

t  time 

t +  transference  number  of  positive  ions 

u  +  effective  mobility  of  positive  ions 

u  ~  effective  mobility  of  negative  ions 

ufr  mobility  of  positive  ions 

ufI“  mobility  of  negative  ions 

v  scalar  convection  velocity  of  positive  and 

negative  ions 
x  position 

Xi  position  of  the  left  end  of  the  porous  electrode  in 

the  circular  cell 

x2  position  of  the  right  end  of  the  porous  electrode 

in  the  circular  cell 

Greek  symbols 

Xt  characteristic  function  of  the  interval  I 

S(  Dirac  delta  functional  centered  at  x  =  £ 

e  porosity  of  the  electrode  material 

<t>  electrostatic  potential 

/tfrec  equivalent  conductance  of  the  electrolyte 

tr  +  effective  ionic  conductivity  of  positive  ions 
cr  “  effective  ionic  conductivity  of  negative  ions 
tr'J*  ionic  conductivity  of  positive  ions 

o-1!?*  ionic  conductivity  of  negative  ions 
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